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Abstract 

We consider nonlinear hyperbolic conservation laws, posed on a dif- 
ferential (n + l)-manifold with boundary referred to as a spacetime, and 
in which the "flux" is deflned as a flux Geld of n-forms depending on a 
parameter (the unknown variable). We introduce a formulation of the ini- 
tial and boundary value problem which is geometric in nature and is more 
natural than the vector field approach recently developed for Riemannian 
manifolds. Our main assumption on the manifold and the flux field is a 
global hyperbolicity condition, which provides a global time-orientation 
as is standard in Lorentzian geometry and general relativity. Assuming 
that the manifold admits a foliation by compact slices, we establish the 
existence of a semi-group of entropy solutions. Moreover, given any two 
hypersurfaces with one lying in the future of the other, we establish a 
"contraction" property which compares two entropy solutions, in a (geo- 
metrically natural) distance equivalent to the distance. To carry out 
the proofs, we rely on a new version of the Unite volume method, which 
only requires the knowledge of the given n-volume form structure on the 
(n -|- l)-manifold and involves the total flux across faces of the elements 
of the triangulations, only, rather than the product of a numerical flux 
times the measure of that face. 

1 Introduction 

The development of the mathematical theory (existence, uniqueness, qualitative 
behavior, approximation) of shock wave solutions to scalar conservation laws 
defined on manifolds is motivated by similar questions arising in compressible 
fluid dynamics. For instance, the shallow water equations of geophysical fluid 
dynamics (for which the background manifold is the Earth or, more generally. 
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Riemannian manifold) and the Einstein-Euler equations in general relativity 
(for which the manifold metric is also part of the unknowns) provide important 
examples where the partial differential equations of interest are naturally posed 
on a (curved) manifold. Scalar conservation laws yield a drastically simplified, 
yet very challenging, mathematical model for understanding nonlinear aspects 
of shock wave propagation on manifolds. 

In the present paper, given a (smooth) differential (n+ l)-manifold M which 
we refer to as a spacetime, we consider the following class of nonlinear conser- 
vation laws 

d{u;{u)) = 0, u = u{x), X G M. (1.1) 

Here, for each m S M, w = uj{u) is a (smooth) Held of n- forms on M which we 
refer to as the flux field of the conservation law (jl.ip . 

Two special cases of (|l.ip were recently studied in the literature. When 
M = R_|_ X N and the n-manifold N is endowed with a Riemannian metric h, 
the conservation law (jl.ip is here equivalent to 

dtu + A\vu{b{u)) =0, u = u{t, y), t>0,ye N. 

Here, divh denotes the divergence operator associated with the metric h. In this 
case, the flux field is a iJux vector field b = b{u) on the n-manifold N and does 
not depend on the time variable. More generally, we may suppose that M is 
endowed with a Lorentzian metric g and, then, (|1.1|) takes the equivalent form 

<l\Vg{a{u)) =0, u ~ u{x), X £ M. 

Observe that the flux a — a{u) is now a vector field on the [n + l)-manifoid M . 

Recall that, in the Riemannian or Lorentzian settings, the theory of weak 
solutions on manifolds was initiated by Ben-Artzi and LeFloch [3] and further 
developed in the follow-up papers by LeFloch and his collaborators [TJ [21 US [121 
[53] ■ Hyperbolic equations on manifolds were also studied by Panov in [53] with 
a vector field standpoint. The actual implementation of a finite volume scheme 
on the sphere was recently realized by Ben-Artzi, Falcovitz, and LeFloch [5]. 

In the present paper, we propose a new approach in which the conservation 
law is written in the form p.ip . that is, the flux uj = a;(u) is defined as a field of 
differential forms of degree n. Hence, no geometric structure is a priori assumed 
on M, and the sole knowledge of the flux field structure is required. The fact 
that the equation (jl.ip is a "conservation law" for the unknown quantity u can 
be understood by expressing Stokes theorem: for sufficiently smooth solutions 
M, at least, the conservation law (jl.ip is equivalent to saying that the total flux 

I uj{u) = 0, UdM, (1.2) 

JdU 

vanishes for every open subset lA with smooth boundary. By relying on the 
conservation law (|l.ip rather than the equivalent expressions in the special cases 
of Riemannian or Lorentzian manifolds, we are able to develop here a theory of 
entropy solutions to conservation laws posed on manifolds, which is technically 
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and conceptually simpler but also provides a significant generalization of earlier 
works. 

Recall that weak solutions to conservation laws contain shock waves and, 
for the sake of uniqueness, the class of such solutions must be restricted by an 
entropy condition (Lax [lOj)- This theory of conservation laws on manifolds 
is a generalization of fundamental works by Kruzkov [TS] , Kuznetsov [TH] , and 
DiPerna [12] who treated equations posed on the (flat) Euclidian space M". 

Our main result in the present paper is a generalization of the formulation 
and convergence of the finite volume method for general conservation law p.ip . 
In turn, we will establish the existence of a semi-group of entropy solutions 
which is contracting in a suitable distance. 

The first difficulty is formulating the initial and boundary problem for (|l.ip 
in the sense of distributions. A weak formulation of the boundary condition 
is proposed which takes into account the nonlinearity and hyperbolicity of the 
equation under consideration. We emphasize that our weak formulation applies 
to an arbitrary differential manifold. However, to proceed with the development 
of the well-posedness theory we then need to impose that the manifold satisfies 
a global hyperbolicity condition, which provides a global time-orientation and 
allow us to distinguish between "future" and "past" directions in the time- 
evolution. This assumption is standard in Lorentzian geometry for applications 
to general relativity. For simplicity in this paper, we then restrict attention to 
the case that the manifold is foliated by compact slices. 

Second, we introduce a new version of the finite volume method (based 
on monotone numerical flux terms). The proposed scheme provides a natural 
discretization of the conservation law (jl.ip , which solely uses the n- volume form 
structure associated with the prescribed flux field uj. 

Third, we derive several stability estimates satisfied by the proposed scheme, 
especially discrete versions of the entropy inequalities. As a corollary, we obtain 
a uniform control of the entropy dissipation measure associated with the scheme, 
which, however, is not sufficient by itself to the compactness of the sequence of 
approximate solutions. 

The above stability estimates are sufficient to show that the sequence of 
approximate solutions generated by the finite volume scheme converges to an 
entropy measure-valued solution in the sense of DiPerna. To conclude our proof, 
we rely on a generalization of DiPerna's uniqueness theorem |12j and conclude 
with the existence of entropy solutions to the corresponding initial value prob- 
lem. 

In the course of this analysis, we also establish a contraction property for 
any two entropy solutions w, v, that is, given two hypersurfaces H, H' such that 
H' lies in the future of H, 



Here, for all reals u, v, the n-form field f2(?I, v) is determined from the given flux 
field uj{u) and can be seen as a generalization (to the spacetime setting) of the 
notion of Kruzkov entropy |u — tJ|. 




(1.3) 
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Recall that DiPerna's measure-valued solutions were used to establish the 
convergence of schemes by Szepessy [26l [28] , Coquel and LeFloch [9l [TOl [11] , 
and Cockburn, Coquel, and LeFloch [HI U\- For many related results and a 
review about the convergence techniques for hyperbolic problems, we refer to 
Tadmor and Tadmor, Rascle, and Bagneiri [3^. Further hyperbolic models 
including also a coupling with elliptic equations and many applications were 
successfully investigated in the works by Kroner [16j . and Eymard, Gallouet, 
and Herbin [13] . For higher-order schemes, see the paper by Kroner, Noelle, and 
Rokyta [T7|. Also, an alternative approach to the convergence of finite volume 
schemes was later proposed by Westdickenberg and Noelle [3T]. Finally, note 
that Kuznetsov's error estimate [S] [5] were recently extended to conservation 
laws on manifolds by LeFloch, Neves, and Okutmustur [22] , 

An outline of the paper is as follows. In Section [2l we introduce our defini- 
tion of entropy solution which includes both initial-boundary data and entropy 
inequalities. The finite volume method is presented in Section [3l and discrete 
stability properties are then established in Section H] The main statements are 
given at the beginning of Section O together with the final step of the conver- 
gence proof. 

2 Conservation laws posed on a spacetime 
2.1 A notion of weak solution 

In this section we assume that M is an oriented, compact, differentiable (n+ 1)- 
manifold with boundary. Given an {n+ l)-form a, its modulus is defined as the 
{n + l)-form 

|a| := \a\dx° A--- Adx", 

where a = adx^ A • ■ • A dx" is written in an oriented frame determined from 
local coordinates x = (x") = (x", . . . , x"). If If is a hypersurface, we denote by 
i = iH '■ H ^ M the canonical injection map, and by i* = is the pull-back 
operator acting on differential forms defined on M . 

On this manifold, we introduce a class of nonlinear hyperbolic equations, as 
follows. 

Definition 2.1. 1. A flux field uj on the (n+ l)-manifold M is a parametrized 
family uj{u) G A"(M) of smooth fields of differential forms of degree n, that 
depends smoothly upon the real parameter u. 

2. The conservation law associated with a flux field uj and with unknown 
■u : M -> R zs 

d{uj{u)) = 0, (2.1) 

where d denotes the exterior derivative operator and, therefore, d(uj{u)^ is a 
field of differential forms of degree (n + 1) on M . 

3. A flux field LU is said to grow at most linearly if for every 1-form p on M 

sup / \p AduU)iu)\ < oo. (2.2) 

ueRJM 
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With the above notation, by introducing local coordinates x = {x") we can 
write for all u G M 

Uj{u) = UJ°'{u) {dx)a, 

{dx)a := A . . . A da;°-i A dx°'+^ A ... A dx". 

Here, the coefficients = lo°'{u) are smooth functions defined in the chosen 
local chart. Recall that the operator d acts on differential forms with arbitrary 
degree and that, given a p-form p and a p'-form /?', one has d{dp) = and 
d{p A p') = dpAp' + {-l)Pp A dp'. 

As it stands, the equation (j2.f p makes sense for unknown functions that are, 
for instance, Lipschitz continuous. However, it is well-known that solutions to 
nonlinear hyperbolic equations need not be continuous and, consequently, we 
need to recast ()2.ip in a weak form. 

Given a smooth solution u of ()2.1|) we can apply Stokes theorem on any open 
subset U that is compactly included in M and has smooth boundary dU. We 
obtain 



= / d{uj{u)) = / i*iuj{u)). (2.3) 

Ju JdU 

Similarly, given any smooth function t/i : Af — > R we can write 

d{ipuj{u)) — dip A Lu{u) + %jj d{uj{u)), 
where the differential dtp is a 1-form field. Provided u satisfies ()2.ip . we find 



d{'ipi Lj{u)) = / dtp /\u){u) 
M Jm 



and, by Stokes theorem. 



dtpAuiu) = / i*{tpuj{u)). (2.4) 

M JdM 

Note that a suitable orientation of the boundary dM is required for this formula 
to hold. This identity is satisfied by every smooth solution to ()2.f |1 and this 
motivates us to reformulate (j2.f p in the following weak form. 

Definition 2.2 (Weak solutions on a spacetime). Given a flux field with at 
most linear grovoth u), a function u g L^{M) is called a weak solution to the 
conservation law \2.1\ posed on the spacetime M if 

dtp A w(?i) = 

M 

for every function i/) : M ^ M compactly supported in the interior M . 

The above definition makes sense since the function u is integrable and Loiu) 
has at most linear growth in u, so that the (n + f )-form dtp A lj{u) is integrable 
on the compact manifold M . 
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2.2 Entropy inequalities 

As is standard for nonlinear hyperbolic problems, weak solution must be further 
constrained by imposing initial, boundary, as well as entropy conditions. 

Definition 2.3. A (smooth) field of n-forms = Q{u) is called a (convex) en- 
tropy flux field for the conservation law p. II) if there exists a (convex) function 
U ■.R->R such that 

n{u) ^ [ duU{v)duUj{v)dv, ueR. 
Jo 

It is said to be also admissible if, moreover, sup \duU\ < oo. 

For instance, if one chooses the function U(u,v) |lt — U|, where w is a real 
parameter, the entropy flux field reads 

fl{u,v) :— sgn{u ~v) {lj{u) — Lj{v)), (2.5) 

which is a generalization to a spacetime of the so-called Kruzkov's entropy- 
entropy flux pair. 

Based on the notion of entropy flux above, we can derive entropy inequalities 
in the following way. Given any smooth solution u to (|2.ip . by multiplying (|2.ip 
by duU{u) we obtain the additional conservation law 

d{n{u)) - {dn){u) + duU{u){dLj){u) = 0. 

However, for discontinuous solutions this identity can not be satisfied as an 
equality and, instead, we should impose that the entropy inequalities 

d{n{u)) ~ {dn){u) + duU{u){duj){u) < (2.6) 

hold in the sense of distributions for all admissible entropy pair {U,il). These 
inequalities can be justified, for instance, via the vanishing viscosity method, 
that is by searching for weak solutions that are realizable as limits of smooth 
solutions to the parabolic regularization of (|2.ip . 

It remains to prescribe initial and boundary conditions. We emphasize that, 
without further assumption on the flux field (to be imposed shortly below), 
points along the boundary dM can not be distinguished and it is natural to 
prescribe the trace of the solution along the whole of the boundary dAI. This 
is possible provided the boundary data, ub ■ dM —^ R, is assumed by the 
solution in a suitably weak sense. Following Dubois and LeFloch [13], we use 
the notation 

e SuM^b) (2.7) 
for all convex entropy pair ([/, 51), where for all reals u 

£u,n{u) := |iJ e M I E{u,v) := n(u) + duU{u){uj{v) - u{u)) < n{v)]. 
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Recall that the boundary conditions for hyperbolic conservation laws (posed on 
the Euclidian space) were first studied by Bardos, Leroux, and Nedelec [3] in the 
class of solutions with bounded variation and, then, in the class of measured- 
valued solutions by Szepessy [57]. Later, a different approach was introduced 
by Cockburn, Coquel, and LeFloch [5] (see, in particular, the discussion p. 701 
therein) in the course of their analysis of the finite volume methods, which was 
later expanded in Kondo and LeFloch [15]. An alternative and also powerful 
approach to the boundary conditions for conservation laws was independently 
introduced by Otto [H] and developed by followers. In the present paper, our 
proposed formulation of the initial and boundary value problem is a generaliza- 
tion of the works [S] and [T5] . 

Definition 2.4 (Entropy solutions on a spacetime with boundary). Let oj = 

ijj{u) he a flux field with at most linear growth and ub G L^{dM) be a prescribed 
boundary function. A function u e L^(Af) is called an entropy solution to the 
boundary value ■problem (|2.ip and \2.7\ if there exists a bounded and measurable 
field of n- forms 7 S L^A"(i9M) such that 

/ (d^j A n{u) + iP {dn){u) - ijduU(u)(duj){u)) 

+ / 'ip\dM{i*^{uB) + duU{uB){l -i*uj{uB))) >0 

JdM 

for every admissible convex entropy pair [U, SI) and every smooth function : 
M R+. 

Observe that the above definition makes sense since each of the terms dV' A 
^(m), ((if2)(u), {dix)){u) belong to L^{M). The above definition can be gener- 
alized to encompass solutions within the much larger class of measure-valued 
mappings. Following DiPerna [l2j, we consider solutions that are no longer 
functions but Young measures, i.e, weakly measurable maps v : M ^ Prob(]R) 
taking values within is the set of probability measures Prob(R). For simplicity, 
we assume that the support supp v is a compact subset of R. 

Definition 2.5. Given a flux field lu = toiu) with at most linear growth and 
given a boundary function ub G L°°{dM), one says that a compactly supported 
Young measure u : M Prob{R) is an entropy measure-valued solution to the 
boundary value problem I \2. 7| ) if there exists a bounded and measurable field 

of n-forms 7 £ L°°A"(9M) such that the inequalities 

J^^ (v, A 17(-) -f i: {d{n{-)) ~ duU{-)iduj){-))) 

^ Id {''^{^*^^'^b) + duU{uB){-i - i*uj{uB)))'^ >0 

hold for all convex entropy pair [U, fJ) and all smooth functions V' ^ 0. 



7 



2.3 Global hyperbolicity and geometric compatibility 

In general relativity, it is a standard assumption that the spacetime should be 
globally hyperbolic. This notion must be adapted to the present setting, since 
we do not have a Lorentzian structure, but solely the n-volume form structure 
associated with the flux field w. 

We assume here that the manifold M is foliated by hypersurfaces, say 

M= U Ht, (2.8) 

0<t<T 

where each slice has the topology of a (smooth) n-manifold N with bound- 
ary. Topologically we have M = [0, T] x N, and the boundary of M can be 
decomposed as 

dM ^ HqUHtUB, 
B = {0,T)xN:^ y dHt. ^^'^^ 

Q<t<T 

The following definition imposes a non-degeneracy condition on the averaged 
flux on the hypersurfaces of the foliation. 

Definition 2.6. Consider a manifold M with a foliation (|2.8ll - (|2.9p and let 

Lo ~ w(?l) be a flux field. Then, the conservation law (|2.ip on the manifold 
M is said to satisfy the global hyperbolicity condition if there exist constants 
< c < c such that for every non-empty hypersurface e G Ht, the integral 
i*(9,jCi;(0) is positive and the function <y9e : R — s- R, 

satisfies 

c<du^e{u)<c, ueR. (2.10) 

The function ipe represents the averaged Hux along the hypersurface e. From 
now we assume that the conditions in Definition 12.61 are satisfied. It is natural 
to refer to Hq as an initial hypersurface and to prescribe an "initial data" 
uq : Hq — > R on this hypersurface and, on the other hand, to impose a boundary 
data ub along the submanifold B. It will be convenient here to use the standard 
terminology of general relativity and to refer to Ht as spacelike hypersurfaces. 

Under the global hyperbolicity condition ()2.8|) - ()2.10p . the initial and bound- 
ary value problem now takes the following form. The boundary condition (|2.7p 
decomposes into an initial data 

uho = "0 (2.11) 

and a boundary condition 

u\g G £u,n{uB). (2.12) 
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Correspondingly, the condition in Definition 12.41 now reads 

^d^ A n{u) + ip {dn){u) ~ipduU{u){duj){u) 



M 

+ / ip\dM {■i*^{uB) + duU{uB){j - i*L0{uB))) 

J B 

+ I l*n{uH^) - I i*Vt{uo) > 0. 

Finally, we introduce: 

Definition 2.7. A flux field lo is called geometry-compatible if it is closed for 
each value of the parameter, 

{duj){u)=Q, ueM.. (2.13) 

This compatibility condition is natural since it ensures that constants are 
trivial solutions to the conservation law, a property shared by many models 
of fluid dynamics (such as the shallow water equations on a curved manifold). 
When (|2.13p holds, then it follows from Definition 12.31 that every entropy flux 
field fl also satisfies the condition 

{dn){u)^0, ueR. 



In turn, the entropy inequalities (j2.6[) for a solution u : AI ^ R simplify drasti- 
cally and take the form 

d{n{u)) < 0. (2.14) 



3 Finite volume method on a spacetime 
3.1 Assumptions and formulation 

From now on we assume that the manifold M — [0, T] x is foliated by slices 
with compact topology N, and the initial data uq is taken to be a bounded 
function. We also assume that the global hyperbolicity condition holds and 
that the flux field uj is geometry-compatible, which simplifies the presentation 
but is not an essential assumption. 

Let T'' = UifeT'' K he a. triangulation of the manifold M, that is, a collec- 
tion of finitely many cells (or elements), determined as the images of polyhedra 
of R"+^, satisfying the following conditions: 

• The boundary dK of an element K is a piecewise smooth, n-manifold, 
dK = Ueca-fs" ^ contains exactly two spacelike faces, denoted by 
and e^, and "vertical" elements 

e°ed°K:=dK\{e+,e],}. 
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• The intersection K C] K' of two distinct elements K, K' 6 is either 
a common face of K, K' or else a submanifold with dimension at most 
(n-1). 

• The triangulation is compatible with the foliation (|2.8p - (|2.9p in the sense 
that there exists a sequence of times ti^ — Q < ti < . . . < tj^j — T such that 
all spacelike faces are submanifolds of := Ht^ for some n — Q, . . . , N ^ 
and determine a triangulation of the slices. We denote by the set of 
all elements K which admit one face belonging to the initial hypersurface 
Ho- 

We define the measure |e| of a hypersurface e C M by 

i*duOj{Q). (3.1) 

This quantity is positive if e is sufficiently "close" to one of the hypersurfaces 
along which we have assumed the hyperbolicity condition (j2.10l) . Provided 
|e| > which is the case if e is included in one of the slices of the foliation, we 
associate to e the function (/jg : R M, as defined earlier. Recall the following 
hyperbolicity condition which holds along the triangulation since the spacelike 
elements are included in the spacelike slices: 

c < du^^± (u) < c, K eT^. (3.2) 

We introduce the finite volume method by formally averaging the conser- 
vation law (|2.ip over each element K € of the triangulation, as follows. 
Applying Stokes theorem with a smooth solution u to (|2.ip . we get 

= / d{u{u)) = I i*u(u). 

JK JdK 

Then, decomposing the boundary dK into its parts e^, e]~-, and ff^K we find 
/ i*Lu{u) - / i*u}{u) + Yj I = 0- (^■^) 

Given the averaged values along e]^ and ^ along e'^ G d'^K, we need an 

approximation uj^ of the average value of the solution u along e^. To this end, 
the second term in (|3.3p can be approximated by 



i*uj{u) « / i*uj{uj^) = |e^.|(^g^(u^) 
and the last term by 

i*Uj{u) « qK,e'>{uK''^K q)' 
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where the total discrete flux qK,e° : ^ R (i-e., a scalar- valued function) must 
be prescribed. 

Finally, the proposed version of the finite volume metliod for the conservation 
law (|2.ip takes the form 

li*u{u+)^li*uj{u]^)- qK,eo{u]^,u]^^J (3.4) 

or, equivalently, 

le^l'^e+l^i) = |eA'l'^e^("K) - lK.eo{uK,u]^^J. (3.5) 

We assume that the functions gx.e" satisfy the following natural assumptions 
for all u, iJ e M : 

• Consistency property : 

qK,eoiu,u)= / i*uj{u). (3.6) 



• Conservation property : 

qK,eoiv,u) ^ -qK^o,e°{u,v). (3.7) 

• Monotonicity property : 

duqK.e'>{u,v) >0, dwqK.,e"{u,v) <0. (3.8) 

We note that, in our notation, there is some ambiguity with the orientation 
of the faces of the triangulation. To complete the definition of the scheme we 
need to specify the discretization of the initial data and we define constant initial 
values uk,o = (for K G Tq) associated with the initial slice Hq by setting 

i*uj{u~^) := i*uj{uo), ej^ C Hq. (3.9) 

Finally, we define a piecewise constant function m'' : A/ ^ R by setting for every 
element K eT^ 

u''{x)^u]^, x€K. (3.10) 

It will be convenient to introduce Nx ^d'^K, the total number of "ver- 
tical" neighbors of an element K G T'', which we suppose to be uniformly 
bounded. For definiteness, we fix a finite family of local charts covering the 
manifold M, and we assume that the parameter h coincides with the largest 
diameter of faces of elements K G T'', where the diameter is computed with 
the Euclidian metric expressed in the chosen local coordinates (which are fixed 
once for all and, of course, overlap in certain regions of the manifold). 
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For the sake of stability we will need to restrict the time-evolution and impose 
the following version of the Courant-Friedrich-Levy condition: for all K G T'', 



Nk 

-T— max sup 

"K 



duuj{u) <infa„(/9g+, (3-11) 



in which the supremum and infinium in u are taken over the range of the initial 
data. 

We then assume the following conditions on the family of triangulations: 



lim ^SS^lZ — = lini ^ = (3.12) 

where Tmax '■— maxi(ii+i — ti) and Tmin '■— mini(<i+i — ti). For instance, these 
conditions are satisfied if Tmax, Tmin, and h vanish at the same order. 

Our main objective in the rest of this paper is to prove the convergence of the 
above scheme towards an entropy solution in the sense defined in the previous 
section. 

3.2 A convex decomposition 

Our analysis of the finite volume method relies on a decomposition of (j3.5p into 
essentially one-dimensional schemes. This technique goes back to Tadmor pS] , 
Coquel and LeFloch [S], and Cockburn, Coquel, and LeFloch [S]. 

By applying Stokes theorem to (|2.13p with an arbitrary m S K, we have 

= / d{uj{u)) = / i*uj{u) 

Jk JdK 

= f^e.m- f .M^)+ E '^-"(^^^)- 

Choosing u ~ u^, we deduce the identity 

|eilV'e+ = |eKlV'e^("K) - <lK.e«{uK,UK), (3.13) 



which can be combined with (|3.5p so that 



By introducing the intermediate values given by 

Nk 

a + 

-K 
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we arrive at the desired convex decomposition 

^4^"^^" 77^ ^ '^ej("K,eo)- (3-15) 

Given any entropy pair {U, fl) and any hypersurface e G M satisfying |e| > 
we introduce the averaged entropy flux along e defined by 

ip^iu) ■.= jj*n{u). 

Obviously, we have ^p'^{u) = ife{u). 

Lemma 3.1. For every convex entropy flux f7 one has 



("i)^^ E ^eV^Lo)- (3-16) 



The proof below will actually show that the function (^^+ o {ip'^+ ) ^ is convex. 

Proof. It suffices to show the inequality for the entropy flux themselves, and 
then to average this inequality over e. So, we need to check: 

Namely, we have 

E (^("x) -^(^i,eo))5nC/(u+) + ^ J2 ^K.eO, 

with 

^1 

DK,eO 9„„L/(w^.)(^(*^A',eO -^J,e")) ^^(^i,eo)) - u'^ ^^o) da. 

In the right-hand side of the above identity, the former term vanishes identically 
in view of p.l4p while the latter term is non- negative since U{u) is convex in u 
and 9„w is a positive n-form. □ 

4 Discrete stability estimates 
4.1 Entropy inequalities 

Using the convex decomposition (j3.15p . we can derive a discrete version of the 
entropy inequalities. 



13 



Lemma 4.1 (Entropy inequalities for the faces). For every convex entropy 
pair {U, f2) and all K £ and G d^K , there exists a family of numerical 
entropy flux functions QK,e° : ^ K satisfying the following conditions for all 
u,v e K; 

• Qk^c" consistent with the entropy flux fi; 

QK.e-{u,u) = / 1*^1{U). (4.1) 

• Conservation property: 



lK,e' 



,o{u,v) = -QK^o,eo{v,u). (4.2) 



• Discrete entropy inequality: with the notation introduced earlier, the finite 
volume scheme satisfies 



(4.3) 

Combining Lemma 13.11 with the above lemma immediately implies: 

Lemma 4.2 (Entropy inequalities for the elements). For each K E one has 
the inequality 

|eil(^4K)-^4("K))+ E {Q{u^,u],J~Q{u]„u-^)) <0. (4.4) 

Proof of Lemma \4.1\ Step 1. For u,w G M and e" G ff^K we introduce the 
notation 



Hx.eO {u, V) ip^+ (u) - (qK,eO (U, «) - qK,e« {u, ")) ■ 



Nk 
Observe that 

We claim that Hj^ satisfies the following properties: 

— 7?K,eO (U, V) > 0, Q^HK,eO (u, v) > 0. (4.5) 

The proof of the second property is immediate by the monotonicity property 
p.Sp . whereas, for the first one, we use the CFL condition (|3.11l) together with 
the monotonicity property (|3.8p . From the definition of Hx i,o{u, v), we observe 
that 

HK,eo{u,UK^o) ^ - E °^K,eo)'Pe+{u) + ^ "i^.e" "^6+ ("K^o ) ' 
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where 

1 <lK.eO (u, UK^o ) - iK.e" 

This gives a convex combination of ip^+{u) and (uK^a)- Indeed, by the 
monotonicity property (|3.8p we have '^eO<^d'>K '^K,e° ^ and the CFL condition 
p. lip gives us 



ClK,e° {U, UK^o ) - QK,e° {u, u) 



< 1. 



Step 2. It is sufficient to estabhsh the entropy inequahties for the family 
of Kruzkov's entropies Jl. In connection with this choice, we introduce the 
numerical version of Kruzkov's entropy flux 

Q(w, V, c) :— qK,e'>{u V c, w V c) — qK.e° {u /\ c,v A c), 

where a V& = max (a, 6) and aAfe = min(a, 6). Observe that Q^.e" (u, "i") satisfies 
the first two properties of the lemma with the entropy flux replaced by the 
Kruzkov's family of entropies 17 = defined in (|2.5p . 
First, we observe 



Hj^^o (uW c,v W c) — Hx^e" [u /\ c,v /\ c) 

= </5^+ (W V C) - {qx.eO [UM C,vy C) - qK,eO [WV c,u'V c)) 



e 



- Up^+ [uAc) ~ — p- {qK,e'> {u Ac,v Ac) - qx^e" (w A c, u A c)) j 
= (u, c) - (q(m, V, c) - Q(u, M, c)) , 

|ej| V y 

where we used 

(f + (-U V c) - ip + (m A c) = / j*r2(M, c) = ip^+ {u, c). 

Second, we check that for u = u^, v = ^ and for any c e R 
-ff^^eo V c, u^^^ Vc) - HK,e«{uK Ac,u^^^ A c) > ip'^+{u+^^„,c). (4.7) 
To prove (|4.7p we observe that 

H^.e" (W) *^) V Hx.e" {^1 A) < Hx.e" (u V A, w V A), 
HK,e" ("i ''^) /\ Hj^ pO (A, A) > Hj(^o (u A A, w A A), 
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where Hx e" is monotone in both variables. Since (z) + is monotone, we have 

HK,eO ("k V C, U]^^^^ V c) - HK,eO ("k ^ C, u]^^^^ A c) 
> i?x_eO ("k , "k^o ) - Hj^^^a (c, c) = (w J gO ) - </:'e+ (c) 

= sgn (u+ ) - (^^+ (c)) (2+^0 ) - ip^+ (c)) 

= sgn(M+^o -c)((P^+(m+ ^0 (c)) = (p^Uux.eO'C). 

Combining this identity with (j4.6p (with u = it^, w = m]^ we obtain the 
following inequality for the Kruzkov's entropies 



'^^("K.eO'C) ~ (■"k'C) + ■|-^(^Q(w,u,c) - Q(w,u,c)) < 0. 

As already noticed, this inequality implies a similar inequality for all convex 
entropy flux flelds and this completes the proof. □ 

If is a convex function, then a modulus of convexity for V is any positive 
real f3 < infl^", where the infimum is taken over the range of data under 
consideration. We have seen in the proof of Lemma [3.11 that (p^ o {(p'^)~^ is 
convex for every spacelike hypersurface e and every convex function U (involved 
in the definition of f2). 

Lemma 4.3 (Entropy balance inequality between two hypersurfaces). For K G 

, let (3^+ be a modulus of convexity for the function (p^\ o ((^'^+ ) and set 
(i ~ min^g7-/i /? + . Then, for i < j one has 



=+ ilrr+ 



^^^^ ^'^'^lU.^,) 

e°Ga'>K (4.8) 



where 7^'' is the subset of all elements K satisfying & Ht^ while ^ -j := 

Uj<fc<j ^t- 

We observe that the numerical entropy flux terms no longer appear in (14. 8p . 

Proof. Consider the discrete entropy inequality (|4.3[) . Multiplying by \e^\/NK 
and summing in if G T^, e° S d^K gives 

E #^^ei("L")- E 141^^. ("k) 

e^ea^if (49) 

+ E {QK,eO {Uk, ) - QK,e» {u^, Ur)) < 0. 
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Next, observe that the conservation property (|4.2p gives 



E QkA^k^^k^o)=^- (4-10) 

e 

So becomes 



(4.11) 



Now, if V is a convex function, and if ^; = ajVj is a convex combination 
of Vj , then an elementary resuh on convex functions gives 

V{v) + f E a,\v, -vf<Y^ a,V{v,), 
j j 

where (3 = inf V" , the infimum being taken over aU Vj. We apply this inequality 
with V = (wj) and V = o ((^"+ which is convex. 

Thus, in view of the convex combination (j3.15|) and by multiplying the above 
inequality by |e^| and then summing in X G T'', we obtain 

e°ea°K e°(id°K 
e'>ed°K 

Combining the result with (|4.11|) . we find 

E \-i\^%i-i<)- E 141^3-. E f^i%o-4P 

KeT'' ATST'' A'gT'' ^ 

- E QK,eo{u](,U]^). 

(4.12) 

Using finally the identity 



J K JdK 



17 



we obtain the desired inequality, after further summation over all of the elements 
K within two arbitrary hypersurfaces. □ 



We apply Lemma l4.3l with a specific choice of entropy function U and obtain 
the following uniform estimate. 

Lemma 4.4 (Global entropy dissipation estimate). The following global esti- 
mate of the entropy dissipation holds: 

E W^KeO-^M'^C [ **"M (4.13) 

for some uniform constant C > 0, which only depends upon the flux field and 
the sup-norm of the initial data, and where fl is the n-form entropy flux field 
associated with the quadratic entropy function U{u) — 

Proof. We apply the inequality (|4.8|) with the choice U{u) ~ 



2 Nk 



After summing up in the "vertical" direction and keeping only the contribution 

-h 




of the elements K G 7^'' on the initial hypersurface i?o, we find 



E ^'3|"i,eo-<r<| E I^kI^^^(-K,0). 



e°ed°K 



Finally, we observe that, for some uniform constant C > 0, 
y \e],\^'^-JuK,o) < C I i*n{uo). 



These expressions are essentially norm of the initial data, and the above 
inequality can be checked by fixing a reference volume form on the initial hy- 
persurface Hq and using the discretization (j3.9p of the initial data uq. □ 



4.2 Global form of the discrete entropy inequalities 

We now derive a global version of the (local) entropy inequality (|4.3p . i.e. we 
obtain a discrete version of the entropy inequalities arising in the very definition 
of entropy solutions. 

One additional notation is necessary to handle "vertical face" of the trian- 
gulation: we fix a reference field of non-degenerate n-forms uj on M which will 
be used to measure the "area" of the faces e° G dK^. This is necessary in our 
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convergence proof, only, not in the formulation of the finite volume method. So, 
for every K £ T'' we define 



|e"|5 := / i*u for faces e'^ G d^'K (4.14) 

and the non-degeneracy condition means that \e'^\[i > 0. 

Given a test-function ?/> defined on M and a face e" G 9°i\r of some element, 
we introduce the following averages 

/ "0 'i'*^ 1 

where, for the first time in our analysis, we use the reference n-volume form uj. 

Lemma 4.5 (Global form of the discrete entropy inequalities). Let 6e a 

convex entropy flux field, and let ip he a non-negative test-function supported 
away from the hypersurface t ^ T. Then, the finite volume approximations 
satisfy the global entropy inequality 

KGTi--^^ iferj* ''■^ (4.15) 



with 



<A'^{ij)+B'\tp)-^CHi^), 



^"i^)--- E ^(^aoK-V^eo)(^^.(Sf+,o)-^^+«)), 



Proof. From the discrete entropy inequalities (|4.3p . we obtain 



(4.16) 



Thanks the conservation property (|4.2p . we have 
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and, from the consistency property (|4.ip . 



e°Gd°K e°ed°K 



■K) 



Next, we observe 

e'>ed°K 



U " J, 



where, we used the inequahty p.l6p and the convex combination p.lSp . In view 
of ^ 

the inequahty (|4.16[) becomes 



< 



I I f 



e°ea°K e°ea°K 

(4.17) 

Note that the first term in (|4.17p can be written as 



Y \^K\'^a°K('P^+{uK) - ip^+{ui^) 



Y I ij{t*i^{u+)-i*n{u],})+ Y I {i'aoK-ij){i*^{u+)~t*rt{u],)). 

^ ^rru Jet. T^^^u Jet. 
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We can sum up (with respect to K) the identities 



K JdK 



Jet Je- _nr^,^JeO 



and combine them with the inequahty (j4.17p . Finally, we arrive at the desired 
conclusion by noting that 

□ 



5 Convergence and well-posedness results 

We are now in a position to establish: 

Theorem 5.1 (Convergence of the finite volume method). Under the assump- 
tions made in Section\^ and provided the flux field is geometry- compatible, the 
family of approximate solutions generated by the finite volume scheme con- 
verges (as h 0) to an entropy solution of the initial value problem (12. 1|) . 
I|2TT|1 . 

Our proof of convergence of the finite volume method can be viewed as a 
generalization to spacetimes of the technique introduced by Cockburn, Coquel 
and LeFloch [BIE] for the (flat) Euclidean setting and already extended to Rie- 
mannian manifolds by Amorim, Ben-Artzi, and LeFloch [T] and to Lorentzian 
manifolds by Amorim, LeFloch, and Okutmustur 

We also deduce that: 

Corollary 5.2 (Well-posedness theory on a spacetime). Let M — [0,r] x N 
be a (n -\- 1)- dimensional spacetime foliated by n-dimensional hypersurfaces Ht 
(t S [0,T]J with compact topology N (cf (12. Let co be a geometry-compatible 
flux field on M satisfying the global hyperbolicity condition (I2.10p . An initial 
data Uq being prescribed on Hq, the initial value problem (|2.ip . (|2.1ip admits 
an entropy solution u € L^[M) which, moreover, has well-defined L^ traces 
on any spacelike hypersurface of M . These solutions determines a (Lipschitz 
continuous) contracting semi-group in the sense that the inequality 



/ i*„,fl{uH',VH') < / i*H^{uH,VH) 

Jh' Jh 



(5.1) 



holds for any two hypersurfaces H, H' such that H' lies in the future of H , and 
the initial condition is assumed in the weak sense 



lim / iH^n{u{t),v{t)) = / iH^ft{uo, 

t>0 Ht J Hq 



vo). (5.2) 
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We can also extend a result originally established by DiPerna [12] (for con- 
servation laws posed on the Euclidian space) within the broad class of entropy 
measure- valued solutions. 

Theorem 5.3. Let uj be a geometry- compatible flux field on a spacetime M 
satisfying the global hyperbolicity condition (j2.10p . Then, any entropy measure- 
valued solution V (see Definition \2.5\) to the initial value problem (j2.1[) . (j2.1ip 
reduces to a Dirac mass at each point, more precisely 

V = (5.3) 

where u G L°°{M) is the unique entropy solution to the same problem. 

We omit the details of the proof, since it is a variant of the Riemannian 
proof given in [1] . 

It remains to provide a proof of Theorem 15.11 Recall that a Young measure 
V allows us to determine all weak-* limits of composite functions a{u^) for all 
continuous functions a, as ft- ^ 0, 

a{u^)^{u,a)^ [ a{X)diy{\). (5.4) 

Lemma 5.4 (Entropy inequalities for the Young measure). Let v be a Young 
measure associated with the finite volume approximations . Then, for every 
convex entropy flux field 17 and every non-negative test-function ^ supported 
away from the hypersurface t = T , one has 

iy,dij An{-)) - [ i*n{uo)<0. (5.5) 

M Jhq 

Based on this lemma, we are now in position to complete the proof of The- 
orem [O] Thanks to (|5.5[) . we have for all convex entropy pairs {U,il), 

d{iy,n{-)) < 

in the sense of distributions on M. On the initial hypersurface Hq the (trace 
of the) Young measure v coincides with the Dirac mass Su,,- By Theorem 15.31 
there exists a unique function u G L°°{M) (the entropy solution to the initial- 
value problem under consideration) such that the measure v coincides with the 
Dirac mass Su- Moreover, this property also implies that the approximations 
converge strongly to u, and this concludes the proof of the convergence of the 
finite volume scheme. 

Proof of Lemma \5.4\ The proof is a direct passage to the limit in the inequality 
(|4.15p . by using the property (|5.4p of the Young measure. First of all, we observe 
that, the left-hand side of the inequality (|4.15p converges to the left-hand side 
of (|5.5p . Indeed, since lo is geometry-compatible, the first term of interest 

V / d{^/jn){u]^) = V / dVAO(u^)== / di;An{u''') 
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converges to jj^^iv, dip A On the other hand, the initial contribution 

in which Uq is the initial discretization of the data uq and converges strongly to 
Mo since the maximal diameter h of the element tends to zero. 

It remains to check that the terms on the right-hand side of (|4.15p vanish 
in the limit h ^ 0. We begin with the first term A'^{i/j). Taking the modulus 
of this expression, applying Cauchy-Schwarz inequality, and finally using the 
global entropy dissipation estimate (|4.13p . we obtain 



e°ea°K 



+ 1 \l/2/ \ 1/2 



hence 



' mm } 



Here, $7 is associated with the quadratic entropy and have used the fact that 
\'4'd°K — ^ C (Tmax + h) . Our conditions (|3.12p imply the upper bound for 
A^{i}}) tends to zero with h. 

Next, we rely on the regularity of "0 and and estimate the second term on 
the right-hand side of (|4.15p . By setting 

we obtain 



e°ea°K 

< 2, SUp|-0eO— '01 / i*^{u]^) — Ceoi*l 



e°ea"K 
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hence 

Again, our assumptions imply the upper bound for B^{ip) tends to zero with h. 
Finally, consider the last term in the right-hand side of (|4.15p 

using the modulus defined in the beginning of Section 2. In view of the inequality 
(|3.17p . we obtain 



and it is now clear that C'* (ip) satisfies the same estimate as the one we derived 
for A''{iP). □ 
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